
double integrand(double x, double param){
  return x * sin(x * param);
}

double primal(double x, double param){
  return sin(x * param) / pow(param, 2)
        - x * cos(x * param) / param;
}
